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1 Introduction 

Finding hints for new physics beyond the standard model (SM) has been a major objective 
of particle physics over the past decades. A very promising strategy to detect such effects 
are high precision experimental measurements which are matched by equally precise the¬ 
oretical predictions. An important ingredient for the precision attainable in a theoretical 
calculation is the knowledge of the coupling constants since they enter the quantum loop 
corrections. 

In this article, we investigate the leading-order hadronic contributions for two of these 
couplings, the electromagnetic fine structure constant, «qed, and the SU(2)l coupling 
constant, ot 2 , both related by the weak mixing angle, 9\y- An accurate knowledge of 
these hadronic contributions is mandatory to accomplish sufficiently precise predictions 
for future high-energy colliders [1] or low energy experiments [2], 

However, the hadronic contributions to the running of «qed turn out to be only poorly 
known at the scale of the Z-boson mass. Compared to oqed at zero momentum transfer, 
there is a five orders of magnitude loss of precision when oqed is taken at the Z-scale 
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turning c>!qed(-^§) into one of the least determined input parameters of the standard 
model [3]. 

Phenomenologically, the leading hadronic contribution to the running of oqed origina¬ 
ting from hadronic vacuum polarisation effects, AaQg D , is determined from a dispersion 
relation and experimental e + e _ scattering data for the hadronic cross-sections [1, 3, 4]. 
Although new data has recently become available, the present analysis does not lead to a 
sufficient improvement of the error which would be needed for the requirements of future 
collider experiments [3]. 

In principle, lattice QCD calculations would be an ideal tool to determine the hadronic 
contributions to electroweak observables such as ckqed or a .2 considered here. However, 
presently the precision that can be obtained from such lattice QCD computations is usu¬ 
ally still lower than from the phenomenological analyses. Nevertheless, the steady progress 
which is taking place in lattice QCD calculations promises to make it an expedient alter¬ 
native to the phenomenological results in the future. In fact, as we will demonstrate 
here, even with our present simulations the statistical uncertainty already matches the 
phenomenological error of AajlS’p and A hvp sin 2 Q\y. 

AQ'qf/d h as been investigated on the lattice for two dynamical twisted mass fermions [5] 
Preliminary results incorporating also dynamical strange and charm quarks for one selected 
momentum value have been reported in [6]. Another determination of Aa™ D , following 
the approach suggested in [4] has been performed in [7]. 

Here, we present our results obtained on the Nf = 2 + 1 + 1 ensembles of the European 
twisted mass collaboration [8, 9]. We will include an estimate of the systematic uncer¬ 
tainties originating from the continuum limit and from the extrapolation to the physical 
point for energies ranging from 0 to lOGeV 2 . 

In contrast to Aa^, the hadronic contributions to the running of the weak mixing angle, 
Owi have not been studied on the lattice so far. Such a calculation is important since the 
phenomenological determination at low energies cannot only be based on data but also 
needs some assumptions such as a partial flavour separation of the cross-section data [10]. 
Lattice calculations can, in contrast, provide a first-principle evaluation of the weak mixing 
angle in the low-momentum region, where several measurements exist [11—14]. In addition, 
due to the great potential of such low energy experiments for unveiling the nature of physics 
beyond the SM, there are also newly planned experimental facilities [15-17], see also [18] 
for a discussion on such experiments. Here, we present the first lattice QCD calculation 
of the leading hadronic contribution to the weak mixing angle, A hvp sin 2 #yE- 

2 The fine structure constant oqed 

Radiative corrections lead to charge renormalisation and thus to the running of the fine 
structure constant obtained by summing the one-particle irreducible bubble insertions in 
the photon propagator [10] 

«QED(0 2 ) - , _ Aa “° ED (Q2) ■ 
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Here, ao is the value at vanishing momentum transfer Q 2 = 0, a 0 1 = 137.035999173(35) [19] 
The leading-order hadronic contribution is given by [1] 

A«QE P D(Q 2 ) = - 4 ™on R (Q 2 ) (2.2) 

and is thus proportional to the subtracted vacuum polarisation function 

n R (Q 2 ) = n(Q 2 ) - n(o). (2.3) 


As mentioned in the introduction, this is usually [1, 10, 19] determined by a phenomenolo¬ 
gical approach relying on the once-subtracted dispersion relation [20] which for Euclidean 
momenta Q 2 reads 

< 2 - 4 > 

and experimental cross-section data for 
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Lattice QCD represents an ab-initio alternative for the calculation of n R (Q 2 ), since the 
hadronic vacuum polarisation tensor can be obtained directly in Euclidean space-time 
from the correlator of two electromagnetic vector currents. 


2.1 Lattice calculation 


The strategy for computing and analysing the hadronic vacuum polarisation function is 
the same as in [22, 23]. In particular, we employ the same set of ensembles [8, 9], which is 
presented in table 1. Additionally, we have checked our chiral extrapolations of the light 
quark contribution by comparing the results with those obtained on a Nf = 2 ensemble 
featuring the physical pion mass [24-26]. The parameters of this ensemble are given in 
table 2. As in [22, 23], we use the conserved point-split vector current at source and sink 
and we restrict our considerations to the quark-connected contributions. In this case, the 
total vacuum polarisation function 

n t°t(Q 2 ) = ^n ud (Q 2 ) + ^n s (Q 2 ) + V(Q 2 ) (2.6) 

is obtained by summing the single-flavour contributions which we define without the charge 
factors. 

For each ensemble and each flavour f, we first fit the temporal vector current correlator to 
determine the vector meson masses, rrq, and their couplings, gi. Then we fit the hadronic 
vacuum polarisation function obtained from the current correlator as detailed in [22] to 
the following functional form 

H f (Q 2 ) = (l - 0(Q 2 - QLtch))nL(Q 2 ) + ©(Q 2 - <3match) n h ig h(<3 2 ), (2.7) 


where @{x) is the Heaviside step function. The low-momentum fit function for 0 < Q 2 < 

Qmatch is g iven by 


nf ow (Q 2 ) = 
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Ensemble 

P 

a[fm] 


7np5[MeV] 

L[fm] 

Aconf 

D15.48 

2.10 

0.061 

48 3 x 96 

227 

2.9 

265/155/156 

D30.48 

2.10 

0.061 

48 3 x 96 

318 

2.9 

203/148/148 

D45.32sc 

2.10 

0.061 

32 3 x 64 

387 

1.9 

397/346/346 

B25.32t 

1.95 

0.078 

32 3 x 64 

274 

2.5 

273/179/180 

B35.32 

1.95 

0.078 

32 3 x 64 

319 

2.5 

201/194/194 

B35.48 

1.95 

0.078 

48 3 X 96 

314 

3.7 

235/103/104 

B55.32 

1.95 

0.078 

32 3 x 64 

393 

2.5 

225/125/125 

B75.32 

1.95 

0.078 

32 3 x 64 

456 

2.5 

158/100/100 

B85.24 

1.95 

0.078 

24 3 x 48 

491 

1.9 

192/142/136 

A30.32 

1.90 

0.086 

32 3 x 64 

283 

2.8 

267/158/158 

A40.32 

1.90 

0.086 

32 3 x 64 

323 

2.8 

248/174/174 

A50.32 

1.90 

0.086 

32 3 x 64 

361 

2.8 

216/147/157 


Table 1. Parameters of the Nf = 2 + 1 + 1 flavour gauge field configurations that have been 
analysed in this work. (3 denotes the gauge coupling, a the lattice spacing, (^) 3 x ^ the space- 
time volume and mps is the value of the light pseudoscalar meson mass. The values for nips 
have been determined in [8]. L is the spatial extent of the lattices. The lattice spacings are taken 
from [21]. The last column gives the number of statistically independent gauge configurations used 
to estimate the light/strange/charm contribution to the polarisation function. 


P 

c S w a[fm] 

(-) 3 x^ 

V a / a 

?np5[MeV] 

L [f'rri] 

-^conf 

2.10 

1.57551 0.094 

48 3 X 96 

128 

4.6 

804 


Table 2. Parameters of ensemble featuring Nf = 2 twisted mass fermions at the physical point. 

and the high-momentum piece for Q^ atch < Q 2 < Q 2 nax reads 

B—l C-l 

n hi g h(Q 2 ) = log(Q 2 ) £ b k (Q 2 ) k + ]T Cl (Q 2 ) 1 . (2.9) 

k= 0 «=0 

The number of terms and thus the fit function is characterised by M, N, B, and C. The 
ansatz in Eq. (2.8) consists of three parts: a series of poles at energies rn t and with residual 
g 2 m 2 , i = 1,... , M, an additive constant oo and further polynomial terms <+ ( Q 2 ) 1 for i > 
1. The poles characterised by ( nii,gi ) are identified with the exponential contributions to 
the time-dependent vector-current 2-point correlation function at zero spatial momentum. 
This is the reason why the parameters ( mi,gt ) are obtained from a fit of the latter and 
inserted into the fit of the vacuum polarisation function under preservation of all error 
correlations. The ansatz in Eq. (2.8) is valid for any four-momentum Q 2 , in particular 
for Q = (Qo,Q = 0^ with Q 2 = Q$- Given such a momentum configuration, the above 
identification follows from the Fourier transform of the polarisation tensor. With the 
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limited statistical precision of the 2-point vector correlator, the number of exponentials 
we can resolve in practice is limited to 2. Contributions from states with even larger 
energies are effectively accounted for by the polynomial terms. 

The ansatz in Eq. (2.9) is chosen to provide an adequate parametrisation of the polarisation 
function. This is the only requirement in the high-momentum region Q^ atch < Q 2 < Q ^ ax . 
While in the low-momentum region the extrapolation beyond the lowest non-zero lattice 
momentum to zero momentum is of physical significance as it predicts the curvature of 
the polarisation function in this interval, in the high-momentum region we only need the 
ansatz in Eq. (2.9) to interpolate the available lattice data. 

In the following, we use Qmax = 100 GeV 2 and choose Qmatch = 2 GeV 2 . Varying Q 2 iatch 
by 1 GeV 2 to the left and to the right gives compatible results. We perform extrapolations 
of the subtracted polarisation function in the light quark mass and the lattice spacing 
only for momenta in the interval 0 < Q 2 < 10 GeV 2 since for larger Q 2 perturbative 
calculations are expected to yield more precise determinations. 

Since our four-flavour ensembles feature unphysically large pion masses, an extrapolation 
to the physical point has to be performed. The pion mass dependence of the single-flavour 
contributions can be assessed by looking at the leading vector meson contribution obtained 
in chiral perturbation theory [27, 28] 

n i (Q 2 )=g 2 v 71 f 2 . (2.10) 

Q z + rriy 

The spectral properties of the heavy vector mesons hardly depend on the pion mass and 
also the coupling constant gy of the p-meson has been found to be well-described by 
a linear fit in the squared pion mass, nip S , cf. [5]. However, the p -meson mass, my , 
strongly depends on the value of the light quark masses, taken to be degenerate in our 
calculation, and thus the squared pion mass [22], This is illustrated in Fig. 1 left by 
a model extrapolation which is constrained by requiring the p-meson mass to attain its 
experimental value [29] at the physical point, similarly to the one used in the two-flavour 
case in Ref. [30]. 

Eq. (??) (2.10) implies that a similar non-linear behaviour can be expected for the light- 
quark hadronic vacuum polarisation function as we indeed observe for the lower set of data 
points in Fig. 1 right. From Eq. (2.10) we also see that we can eliminate this non-linear 
dependence on the squared pion mass to a large extent by employing the lattice redefinition 
presented in Ref. [5] for the light-quark contribution to the vacuum polarisation function, 


5jR (<3 2 ) = n R d 



( 2 . 11 ) 


if we use H = my , i.e. the p-meson mass at unphysically large up and down quark masses. 
The beneficial effect this has on the data points is depicted as the upper set of data points 
in Fig. 1 right. Hence, in the following we use the above redefinition in the light-quark 
sector. For the contributions of the heavy quark flavours, we use the standard definition 
of the vacuum polarisation function n^(Q 2 ), f = s, c. Our redefinition of the total 
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Figure 1. Left: Dependence of p-meson mass, my, on the squared pion mass, m| g , determined 
on the Nf = 2 + 1 + 1 twisted mass ensembles shown in table 1. Right: Chiral extrapolation 
of the renornralised light-quark vacuum polarisation function at a generic value of Q 2 = 1 GeV 2 
obtained with the redefinition given in Eq. (2.11) (upper set of data points) and with the standard 
definition Eq. (2.3) (lower set of data points). For the latter the same model extrapolation as for 
the p-meson mass has been used. 

then follows from the sum over all quark flavours as in Eq. (2.6) 

= —47ra 0 IIp d + g n R ( Q 2 ) + g Hr {Q 2 )^ ■ (2.12) 

The lattice data obtained with the definition in Eq. (2.12) can be sufficiently well 

described already by a linear dependence on the squared pion mass. Since we use only 
this definition throughout this work and no confusion is possible, we henceforth omit the 
bar in AaJ D . 

2.2 Results 

In order to show that the above redefinition in Eq. (2.11) indeed provides the expected 
benefit for the chiral extrapolation of the light quark contribution to the running of the 
fine structure constant, we show the data for both Eqs. (2.2) and (2.11) with H = my 
in Fig. 2 left for a single momentum value Q 2 = lGeV 2 . The upper set of data points 
obtained with the redefinition Eq. (2.11) evidently is much easier to extrapolate to the 
physical value of the pion mass than the lower points procured from the standard definition 
Eq. (2.2). 

Since we do neither observe lattice spacing artefacts nor finite size effects in these data at 
Q 2 = 1 GeV 2 , we can actually compare our results computed on the four-flavour ensembles 
linearly extrapolated in the squared pion mass, mp S , with those obtained from the Nf = 2 
ensemble featuring the physical pion mass. Additionally to the standard analysis, we have 
performed a correlated [1,1] Pade fit [31], possessing the same number of parameters, up 
to Qmax = 1-5 GeV 2 such that Q 2 = 1 GeV 2 is safely covered. As expected, the values for 
the pole parameters determined from the temporal correlator in our standard approach 
and from the Pade fit are compatible 

a 2 rriy = 0.153(35) b n = 0.1575(81) (2.13) 
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Figure 2. Left: Light-quark contribution to Aaq V g D with filled symbols representing points ob¬ 
tained with Eq. (2.11) using H = my, open symbols refer to those obtained with Eq. (2.2), 
i. e. H = 1 in Eq. (2.11). In particular, the two-flavour results at the physical point have been 
computed with the standard definition. The light grey errorband displays the uncertainty of the li¬ 
near fit represented by the black dotted line, wheras the dark grey errorband belongs to a quadratic 
fit represented by the green solid line. Right: Nf = 2 + 1 + 1 contribution to for the three 

lattice spacings at a fixed pion mass of Tops « 320 MeV. 


and also the results of both analyses of the leading hadronic contribution to the running 
of the fine structure constant at the physical point completely agree with each other 
and with the extrapolated result obtained on the four-flavour ensembles indicating that 
the systematic uncertainty caused by the chiral extrapolation is small. The results at 
the physical value of the pion mass, which are depicted in the left panel of Fig. 2, are 
summarised in table 3. 


Nf = 2 + 1 + 1 extrapolated 

Nf = 2 standard Nf = 2 [1,1] Pade 

0.003068(50) 

0.003097(88) 0.003062(77) 


Table 3. Comparison of the chirally extrapolated result for Aa^o^l GeV 2 ) obtained on the 
Nf = 2 + 1 + 1 ensembles with those obtained on the Nf = 2 ensemble at the physical point. For 
the latter, we have performed our standard analysis but without the redefinition and also tested a 
[1,1] Pade fit. 

Including also the heavy quark contributions by using Eq. (2.12), a dependence on the 
lattice spacing is clearly visible, especially in the high-Q 2 region shown in the right panel 
of Fig. 2. This is accounted for by combining the chiral extrapolation with taking the con¬ 
tinuum limit and employing the following fit function to the four-flavour results obtained 
on individual ensembles 

Aa QED( ( 5 2 )( m PS: a) = A + B rrip S + C a 2 (2.14) 

with fit parameters A, B, C for each momentum value Q 2 G {0, 0.02, 0.04,..., 10} GeV 2 . 
In [32] we have shown that automatic O(a) improvement is at work for our definition of the 
hadronic vacuum polarisation function. Thus, performing the continuum extrapolation 
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without a term linear in the lattice spacing a in Eq. (2.14) is justified. The ansatz in 
Eq. (2.14) neglects any dependence on the finite lattice extent which has been found to 
be smaller than our current statistical uncertainties and will be discussed when assessing 
the systematic effects of our calculations. 



Figure 3. Nf = 2 + 1 + 1 contribution to Acq^-p compared to the data collected in [33] employing 
the dispersion relation in Eq. (2.4). The lattice data are taken at the physical point and in the 
continuum limit (CL). 

The results are depicted in Fig. 3 together with the results obtained by a phenomenological 
analysis [1]. Here, for both the lattice calculation and the phenomenological analysis only 
the statistical errors are shown. Over the whole momentum range, perfect agreement with 
comparable statistical uncertainties is found. We will discuss the systematic uncertainties 
of our lattice QCD determination below. An updated phenomenological analysis including 
all data published till the end of 2014 will soon be available [34], The lattice data also 
agree with those results featuring even smaller uncertainties. 

2.2.1 Systematic uncertainty from the choice of vector meson fit ranges 

As mentioned before, the first step in our analysis is the determination of the masses 
and the coupling constants of the vector mesons from the vector two-point functions at 
zero momentum. The values of the spectral parameters differ when varying the fit range. 
We have repeated the complete analysis for various vector meson fit ranges for the light, 
strange and charm quark currents propagating the uncertainty to the final results. 

In the light quark sector depicted in Fig. 4, we observe systematic uncertainties depen¬ 
ding on whether we start fitting the vector meson correlator at 0.6 fm or at 0.7 fm whereas 
changing the upper border of the fit interval by 0.1 fm does not lead to observable ef¬ 
fects. The dependence on the lower starting point of the fit can be attributed to excited 
state contamination of the p-meson correlator. When stating the final results for selected 
momentum values below, we take for these systematic uncertainties half the difference 
between the central values that are furthest apart from each other. 










Figure 4. Light quark contribution to obtained with different fit ranges for the p-meson 

mass, my and coupling, gy. The standard fit range is [0.7fm, 1.2 fm]. 


For the heavy flavours changing the fit interval by 0.1 fm to the left and to the right of both 
the lower and the upper time slice of the fit ranges does not lead to observable differences. 
This is shown in Fig. 5. 
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Figure 5. Dependence of the single-flavour contributions to Aoqed on the fit range of the ss- 
correlator (left panel) and of the J/'L-correlator (right panel). The standard ss-correlator fit range 
is [0.9fm, 1.4fm], whereas the one for the charm quark correlator is [1.2fm, 1.7fm]. The minor 
discontinuities at Q 2 = 2 GeV 2 arise from connecting the low-momentum Eq. (2.8) and high- 
momentum Eq. (2.9) fit functions at this point by a simple step function as shown in Eq. (2.7). 
Due to the subdominance of the heavy flavour contributions, those discontinuities do not influence 
the final result. 


2.2.2 Systematic uncertainty from the choice of vacuum polarisation fit func¬ 
tion 

Performing the whole analysis with different numbers of terms in our vacuum polarisation 
fit functions also leads to observable differences in the light quark contribution as shown 
in Fig. 6. These are larger than the effects from the fit ranges of the vector meson fits dis¬ 
cussed in the preceding subsection and thus present the dominant systematic uncertainty 
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in our calculation. It might be possible to improve the situation by e.g. the method of ana¬ 
lytic continuation [35, 36] or by taking momentum derivatives of the vacuum polarisation 
function [37]. 
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Figure 6. Light quark contribution to Aqq[? i;) obtained from different fit functions. The standard 
fit is M1N2B4C1. 


The situation for the heavy quarks is shown in Fig. 7. Here, almost no systematic devia¬ 
tions are visible. Furthermore the contributions from the heavy quarks are about an order 
of magnitude smaller than the light-quark one. Hence, we do not take systematic effects 
from the variation of the second-generation quark fit functions into account in our final 
error estimate. 
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Figure 7. Dependence of the single-flavour contributions to Aqqed on the choice of fit function 
for the strange (left panel) and for the charm (right panel) quark pieces. For the strange quark 
the standard fit is M1N2B4C1, whereas the one for the charm quark correlator is M1N2B3C5. 


2.2.3 Finite size effects 

In lattice QCD, typically raps L > 4 is required to minimise systematic effects due to the 
finite lattice volumes, where L denotes the spatial extent of the lattice. The Nf = 2 + 1 + 1 
ensembles analysed in this work feature 3.35 < raps L < 5.93. Restricting our data to the 
condition rapsT > 3.8 yields the picture shown in the left panel of Fig. 8. Hence, we do 
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not associate a systematic uncertainty to the usage of ensembles possessing smaller raps-L 
values. 
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Figure 8. Four-flavour contribution to Aa^ D obtained with (left panel) mpsL > 3.35 [standard] 
and mpsL > 3.8 and (right panel) tops < 491 MeV [standard] and mpgL < 400 MeV. 


2.2.4 Systematic uncertainty from including heavy pion masses 

In order to extrapolate to the physical point, m n ~ 140 MeV, often not too high pion 
masses should be included in the fit. The ensembles entering the standard analysis com¬ 
prise pion masses up to mps ~ 491 MeV. Using only the ensembles with mps < 400 MeV 
yields fully compatible results for as can be seen in the right panel of Fig. 8. 

Therefore, we do not account for a systematic uncertainty related to the usage of pion 
masses above 400 MeV. 

2.2.5 Final results for selected momentum values 

Table 2.2.5 contains our final results compared to those of a phenomenological analysis [1] 
utilising the once-subtracted dispersion relation Eq. (2.4). The first error denotes the 
statistical and the second error the systematic uncertainty of our results. The latter 
constitutes the dominant source of uncertainty of which the biggest part originates from 
the choice of the vacuum polarisation fit function. This might change when lowering the 
statistical uncertainty, because then the vacuum polarisation fit gets more constrained. 
Alternatively, avoiding to fit the vacuum polarisation might be considered. 


3 The weak mixing angle sin 2 9\y 


The weak mixing or Weinberg angle, 6\Vi is one of the fundamental parameters of the 
electroweak standard model defined by 


sin 2 9\y = 


9 


/ 2 


g' 2 + g z 


e 2 _ oqed 

g 2 a 2 


(3.1) 


where g is the SJJ{2)l coupling constant and g' the U(l)y coupling constant. The second 
equality is the electroweak unification condition e 2 = g 2 sin 2 9\y for the positron charge e. 
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Q 2 [GeV 2 ] 

this work dispersive analysis [1] 

0.02 

0.163 (05) (09) • 10“ 3 0.174 (02) • 10“ 3 

1.00 

3.721 (96)(145) ■ 10“ 3 3.651 (40)•10“ 3 

2.00 

4.993 (102)(144) •10~ 3 4.916(61) •10“ 3 

3.00 

5.800 (111)(151) • 10“ 3 5.725(74)•10“ 3 

4.00 

6.396 (108)(156)•10“ 3 6.333(84)•10“ 3 

6.00 

7.264 (114) (159) • 10" 3 7.223 (98) • 10“ 3 

8.00 

7.906 (124) (151) • 10“ 3 7.850 (107) • 10" 3 

10.0 

8.419 (130)(159) • 10“ 3 8.420(114)•10“ 3 


Table 4. We tabulate Ao;q'^ > d (Q 2 ) for selected values of Q 2 . The first error of the lattice results is 
statistical, the second systematic. The phenomenological values of Aag'^’j-^Q 2 ) have been obtained 
from the dispersive analysis of Ref. [1]. 

Thus, the running of the weak mixing angle can be obtained from the running of the fine 
structure constant and the SU(2)l coupling a 2 . In the leading logarithmic approximation 
this is given by [38] 

sin 2 6 w (Q 2 ) = sin 2 (fl°) 1 Aa2 (Q ) = s i n 2 ( 6 » 0 )(l + A(Q 2 )) (3.2) 

1 - Ackqed [Q z > 

where sin 2 (0o) = tf, andA(Q 2 ) = Ao!qed(Q 2 ) —AfWQ 2 ) is an abbreviation for A sin 2 0w(Q 2 )- 

a 2 ' ' 

The value of sin 2 ( 6 * 0 ) has essentially been measured by the Boulder group studying atomic 
parity violation in Cesium [11], the latest value is sin 2 (0o) = 0.2356(20) [39]. The standard 
model prediction in the MS scheme is sin 2 (@o) = 0.23871(9) [18, 40] which is the value 
employed in the analysis below in order to gain fully theoretical results without experi¬ 
mental input. In the computation of this value the Higgs boson mass determined by the 
LHC experiments [41, 42] has been used. 

A phenomenological value of the leading hadronic contribution to the running of sin 2 9\y 
between 0 and the Z-scale has been computed for the first time in [43] relying on results 
of [44] . The method has been described and used with an older dispersive analysis [45] be¬ 
fore [46] . In [40] the error has been reduced with respect to the original rather conservative 
estimate of the uncertainty by about an order of magnitude. 

The leading hadronic contribution to the running of the SU(2)l coupling constant origi¬ 
nates from Z -7 mixing 


-had'^/'v/x/v- 

Z 7 
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From the expressions for the hadronic currents of up-type (u) and down-type (d) quarks 


Jl = Jj-sin \e w )j; (3.3) 

= 4 E - 75)«/ - 3/7/i(! - 75)rf/) (3.4) 

/ 

J /I = E - ^5/7/^/) (3-5) 

where 3 refers to the third component of the weak isospin current and 7 to the electro¬ 
magnetic current discussed already above, we see that to leading order 

n Z7 « n 37 = < J 3 J 7 ) (3.6) 

and thus the leading hadronic contribution to the running of 012 is given by [ 1 , 10 ] 

a« 2 VP (Q 2 ) = -g 2 (n 37 (Q 2 ) - n 37 (o)) . (3.7) 


As for the purely electromagnetic current correlator, II 37 denotes the transverse part of 
the vacuum polarisation function. 

Beyond the leading log approximation, AcqwQ and Aa!, lvp become renornralisation scheme 
dependent. Additional hadronic contributions to these corrections at the scale of the W- 
rnass and the Z-mass originate from chiral symmetry breaking. They have been shown to 
be calculable in perturbation theory and to be at least two orders of magnitude smaller 
and thus negligible compared to the leading contributions [10]. Thus, having computed 
A«qed before, all that is left to do to leading order is to compute Act 2 as given in Eq. (3.7). 


3.1 Lattice calculation 


Since our ensembles feature mass-degenerate up and down quarks, m u = rrid, light-quark 
disconnected contributions cannot occur in II 37 due to the isospin symmetry of the vac¬ 
uum. Without those interference terms, single-flavour contributions to Eq. (3.6) in the 
continuum limit have the general structure IT 37 7 =< (V — A)V >, where V and A denote 
vector and axial vector currents, respectively. Since QCD conserves parity, no mixing 
between vector and axial vector currents occurs such that without quark-disconnected 
contributions we obtain for up-type quarks twice the contribution of down-type quarks 

U %’ U = e E<(“/W) («/ 7 ,«/)> = ' ( 3 - 8 ) 


Combining this with AoqeDj the leading-order hadronic contribution to the running of 
the weak mixing angle from the two light flavours reads 


Al7« 2 ) = -A“‘ Q $ vp (<? 2 ) + A" d c^ D (Q 2 ) = J 9 2 n” 


'(Q 2 ) - -e 2 n““{(5 2 ). (3.9) 


Neglecting disconnected contributions also for the heavy flavours, we have for the strange 
and the charm quark contributions 


A s hvp (Q 2 ) = —A s a% vp (Q 2 ) + A s agP D (Q 2 ) = ^ 2 IF S (Q 2 ) - 4 e 2 IT s (Q 2 ) 
Ac Vp (Q 2 ) = -A c 4 vp (Q 2 ) + A c «5 e p d (Q 2 ) = lg 2 n CC (Q 2 ) - §e 2 n~(Q 2 ) , 
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(3.10) 

(3.11) 


0.02 



Figure 9. Nf = 2 + 1 + 1 contribution to Aa^ compared to the data collected in [33] for all 
quarks except the top. The lattice data are extrapolated to the physical point and to the continuum 
limit (CL). 

respectively. Hence, the single flavour contributions are all proportional to the hadronic 
vacuum polarisation function but with different prefactors than for Aa^-p. In order to 
treat both contributions to A hvp sin Q\y consistently, we use in the light sector the same 
redefinition of the vacuum polarisation function for Aa^ vp as for Aa^-p. Thus, in our 
lattice calculation we consider 

A hvp sin 2 0 W {Q 2 ) = A hvp ’ ud (Q 2 • H 2 /H 2 hys ) + A hvp ’ s (Q 2 ) + A hvp ’ c (Q 2 ). (3.12) 

3.2 Results 
3.2.1 A«2 VP 

As stated above, the leading hadronic contribution to the running of the weak mixing 
angle in the leading logarithmic approximation is obtained from the difference of the 
corresponding contributions of the electromagnetic and the SU{2)l coupling constants, 
aQED and « 2 - In contrast to Aaq^-Q, it is not straightforward to extract Aa^' 1 ’ from 
experimental e + e~ —> hadrons data, since the data comprising the three lightest quarks 
would have to be separated either in up-type (u) and down-type (d and s) quarks or 
assuming isospin symmetry in light and strange quark contributions. This problem has 
no unique solution, e.g. final states involving kaons could either originate directly from a 
strange quark current or from a gluon that could be radiated off light quarks. Another 
possibility is to assume SU(3)f symmetry and thus only split the data into information 
attributed to the three lightest quarks and the rest. The contributions from charm and 
heavier quarks can be computed in perturbation theory. 

Fig. 9 shows our results after combined extrapolation to the physical point and to van¬ 
ishing lattice spacing compared to the results of [1], There, two ways of flavour sepa¬ 
ration have been implemented, one is assuming approximate SU (3) f and the other one 
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Nf = 2 + 1 + 1 extrapolated 

Nf = 2 standard Nf = 2 [1,1] Pade 

-0.002717(43) 

-0.002742(78) -0.002710(68) 


Table 5. Comparison of results for A hvp,ud sin 2 9w{l GeV 2 ) at the physical point. The same 
analyses as indicated below table 3 have been performed. 

517(2)/ symmetry neglecting OZI violating terms. Our results clearly prefer the 517(3)/ 
flavour separation and thus indicate that the latter assumption is not tenable as has 
also been observed in [47] in a different context. In Fig. 9 we have multiplied the data 
from [33] with sin 2 9\y(Mz)/ sin 2 (#o) to account for the different reference values em¬ 
ployed. As mentioned before sin 2 (#o) = 0.23871(9) and the value used by Jegerlehner is 
sin 2 9w{Mz) = 0.23153(16) which has been measured at LEP [48]. The flavour separation 
performed for the data set including very recent e + e _ measurements is based on isospin 
symmetry relations [34] and the results are much closer to the ones based on 517(3)/ 
flavour separation in Fig. 9 than to the old 517(2)/ curve. Thus, our lattice results are 
also compatible with the newest phenomenological analysis based on an isospin 517(2)/ 
flavour separation, however, not assuming flavour non-diagonal elements to be small. 

3.2.2 A hvp sin 2 6w 

Having determined the four-flavour contributions to Aa^-p and Aa^, it is straightfor¬ 
ward to obtain the leading-order hadronic vacuum polarisation contribution to the running 
of the weak mixing angle 

A hvp sin 2 6 W (Q 2 ) = Aag£ D (Q 2 ) - AaiT p (Q 2 ) ■ (3-13) 

This is the central observable measured in various low-energy experiments in order to gain 
hints on beyond the SM physics. In subsection 3.2.4 below, a selection of such experiments 
operating at momentum transfers investigated in this work will be listed. 

The physical results for the light-quark contribution for each momentum value can again 
be obtained from extrapolations in the squared pion mass as shown in Fig. 10 for Q 2 = 
1 GeV 2 . In contrast to the case of <uqed depicted in the left panel of Fig. 2, for the weak 
mixing angle combining the redefinitions according to Eq. (2.11) of «qed and «2 leads 
to lower values than obtained with the standard definitions. The common feature of the 
leading-order hadronic contributions of both quantities is that the values procured with 
the redefinitions can be already well-described by a simple linear extrapolation in the 
squared pion mass to the physical point yielding a result which is compatible with those 
of the standard analysis as well as the one from Pade approximants on the ensemble of 
two dynamical quarks at the physical point. The results at the physical value of the pion 
mass are given in table 5. 

When incorporating the heavy quarks, the chiral extrapolation is again combined with 
taking the continuum limit of the four-flavour result according to 

A hvp sin 2 9\y(Q 2 )(mps, a) = A + B mp S + C a 2 . (3-14) 
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mis [GeV 2 ] 

Figure 10. Light-quark contribution to A hvp sin 2 6>w with filled symbols representing points 
obtained with Eq. (2.11), open symbols refer to those obtained with Eqs. (2.2) and (3.7). In 
particular, the two-flavour results at the physical point have been computed with the standard 
definitions. The light grey errorband displays the uncertainty of the linear fit represented by the 
black dotted line wheras the dark grey errorband belongs to the quadratic fit shown as the green 
solid line. 

The results are shown in Fig. 11. Complying with the indication from the previous subsec¬ 
tion, we have employed the results for Aa!, 1 ' 1 ’ obtained from SU (3 )f flavour separation in 
Fig. 11 together with the factor needed to take the different reference values into account. 
Since we do not have information on the correlation of the data in [33], we have simply 
added the uncertainties of Aqq^ and Aa 2 lvp in quadrature and may thus overestimate 
the errors of the phenomenological determination. 

3.2.3 Systematic uncertainties 

Since the systematic uncertainties stem from the same sources as for Ace^-p discussed 
before, the relative errors are the same and only the absolute numbers differ due to the 
different prefactors of the renormalised vacuum polarisation function. Naturally, also the 
plots all look very similar. Therefore, we refrain from discussing the systematic effects 
separately and only summarise the general findings. 

As before, due to the light quark contribution being an order of magnitude bigger than 
the contributions from the heavy quarks, we only need to take systematic uncertainties 
of this part into account. The dominant source of systematic errors is again the choice of 
the vacuum polarisation fit function as depicted in Fig. 12. 

The only other relevant effect comes from the excited state contamination of the p-meson 
correlator and is shown in Fig. 13. Finite volume effects and the choice of rather heavy pion 
masses in the chiral extrapolation seem to be negligible in our calculation as outlined be- 


Nf = 2 result, standard lit ■—B- 

Nf = 2 result, Pade fit >— 
a = 0.086 fm, L = 2.8 fm >—♦- 
a = 0.078 fm, L = 2.5 fm ■— •- 
a = 0.078 fm, L = 1.9 fm — 
a = 0.078 fm, L = 3.7 fm —» 
a = 0.061 fm, L = 1.9 fm —r- 
a = 0.061 fm, L = 2.9 fm 



> 

CD 

0 


<1 
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Figure 11. Nf = 2 + 1 + 1 contribution to the leading-order hadronic contribution A hvp sin 2 9w 
compared to the difference of the data collected in [33]. The lattice data are extrapolated to the 
physical point and to the continuum limit (CL). 
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Figure 12. Light quark contribution to A hvp sin 2 9w obtained from different fit functions. The 
standard fit is M1N2B4C1. 


fore. The only unknown systematic effect is the heavy-flavour disconnected contributions 
which we have neglected here. 

3.2.4 Final results for selected momentum values 

In table 6 we collect our results for A hvp sin 2 6\y with statistical as well as systematic 
uncertainties for selected momentum values. Experiments which have measured or will 
measure the weak mixing angle in the respective momentum region are also indicated. 
The outcome of the E158 experiment at the SLAC linear accelerator was the first suc¬ 
cessful measurement of parity violation in electron-electron (Mpller) scattering [12]. The 
momentum transfer was Q 2 = 0.026 GeV 2 . The Qweak experiment conducted at JLAB in 
2012 measured parity violation in electron-proton scattering at almost exactly the same 
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Figure 13. Light quark contribution to A hvp sin 2 9w obtained with different fit ranges for the p 
meson properties. The standard fit range is [0.7 frn, 1.2 fm]. 


Q 2 [GeV 2 ] 

this work 

experiment 

0.02 

-0.158 (05) (08) • 10" 3 

E158, Qweak 

1.00 

-3.706 (83) (127) • 10“ 3 

PVDIS 

2.00 

-5.021 (96) (135) • 10“ 3 

PVDIS 

3.00 

-5.801 (104) (135) • 10" 3 

Solid 

4.00 

-6.398(102) (135) • 10~ 3 

Solid 

6.00 

-7.251 (111) (136) • 10” 3 

Solid 

8.00 

-7.867(112) (137) • 10~ 3 

Solid 

10.0 

-8.352(119) (138) • 10" 3 

Solid 


Table 6. We tabulate A hvp sin 2 9w for selected values of Q 2 . The first error of the lattice results 
is statistical, the second systematic. Several low energy experiments sensitive to the respective 
momentum regions are indicated in the last column. 


momentum transfer [14], The data is still being analysed. The predicted final uncer¬ 
tainty is about 5% or 0.7 • 10 -3 taking the central SM value. Another JLAB experiment 
performed by the PVDIS collaboration determined the weak mixing angle from parity- 
violating deep inelastic scattering [13, 49] which effectively means electron-quark scatter¬ 
ing at Q 2 = 1.085 GeV 2 and Q 2 = 1.901 GeV 2 . The envisioned successor of the PVDIS 
experiment which also measures parity violation in electron-quark scattering is the SoLID 
spectrometer proposed at JLAB [16]. It can study about 20 kinematic points with Q 2 
ranging from about 2GeV 2 to about 10GeV 2 . Its target accuracy is 6 • 10~ 4 . Our re¬ 
sults in table 6 indicate that it will be essential to take the leading QCD corrections into 
account in order to deploy the whole potential of the experiment in the search for new 
physics beyond the SM. 
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4 Summary and Outlook 


Hadronic contributions to the running of electroweak parameters nowadays constitute the 
major uncertainties of their values even at high energies thus also limiting the precision 
achievable in predictions for future high-energy colliders. Here we have considered the 
running of «qed and of the weak mixing angle which represents one of the most important 
parameters of the SM and provides a sensitive probe of new physics over a large energy 
range. 

Lattice QCD provides a most valuable tool to compute these hadronic contributions from 
first principles alone. As we have demonstrated in this article, lattice QCD can be used 
to compute to a good precision the leading-order hadronic contribution to the running of 
a QED- In particular, we have carried out the first dynamical four-flavour calculation of 
the leading-order hadronic contribution to the running of the fine structure constant and 
the first lattice QCD calculation of the leading hadronic contribution to the shift of the 
weak mixing angle at energies between 0 and lOGeV 2 . In both cases the chiral as well 
as continuum extrapolations have been performed. A main effort has been undertaken to 
assess systematic uncertainties on a quantitative level. 

For both quantities, agreement of our results with a phenomenological determination is 
observed with an even comparable statistical uncertainty. However, we have found that 
the systematic effects of the calculation still exceed the statistical errors. The dominant 
systematic uncertainty has been found to be the choice of fit function. Thus, methods 
which try to avoid fitting the vacuum polarisation are promising to reduce the overall 
uncertainty. A further improvement can be achieved by increasing the statistical precision 
which would more strongly constrain the vacuum polarisation fit. Such improvements 
might be accomplished by the use of the all-mode-averaging [50] or the exact deflation [51, 
52] techniques. 

In this article, we have provided further successful examples for the programme to de¬ 
termine hadronic contributions to electroweak observables from lattice QCD. The steady 
progress in lattice QCD with ever increasing statistical accuracy and better understand¬ 
ing and control of systematic uncertainties makes the lattice approach to compute these 
hadronic contributions very promising and gives hope that lattice results can directly be 
used for future low energy and collider experiments. 
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